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1 Introduction 



A theoretically clean measurement of the angle 7 = arg 




(also denoted as 



03 in the literature) can be obtained using CP- violating B — )■ D^*^K^*^ decays. The 
interference between the h — )■ cus and h — )■ ucs tree amplitudes results in observables 



the relative strong phase 5b between the two amplitudes. The hadronic parameters, 
tb and 5^, depend on the B decay under investigation; they can not be precisely 
calculated from theory, but can be extracted directly from data by simultaneously 
reconstructing several different D final states. 

In this contribution we present the most recent 7 determinations obtained by 
BABAR, based on the full sample (~ 468 x 10^ B"^ decays) of charged B mesons pro- 
duced in e+e- T{4.S) B+B' and accumulated in the years 1999-2007. The fol- 
lowing decays have been reconstructed: (i) B^ — )> D^*'^K^ and B^ DK*^{K*^ — )■ 
K^s^^), with D K^^h+h-, h = tt,K; (ii) B^ DK^, with D decaying to CP- 
eigenstates fcp', (hi) -B^ — )■ D^*^K^, with D decaying to K'^n^. The results are 
statistically limited, as the effects that are being searched for are tiny, since: (i) 
the branching fractions of the B meson decays considered here are on the order of 
5 X 10""^ or lower; (ii) the branching fractions for D^*^ decays, including secondary 
decays, range between O(10~^) and O(10~*^); (iii) the interference between the b ^ c 
and b u mediated B decay amplitudes is low, as the ratios tb are around 0.1 due 
to CKM factors and the additional color-suppression of A{b u). 

The B decay final states are completely reconstructed, with efficiencies between 
40% (for low-multiplicity, low-background decay modes) and 5% (for high-multiplicity 
decays). The selection is optimized to maximise the statistical sensitivity S/\^S + B, 
where the number of expected signal (5*) and background (B) events is estimated 
from simulated samples and data control samples. Signal B decays are distinguished 
from BB and continuum qq background by means of maximum likelihood fits to 
two variables exploiting the kinematic constraint from the known beam energies: 



the energy-substituted invariant mass mEs=y E^ls^m ~ P*b the energy difference 

AE=E*Q—E^^^^^^^. Additional continuum background discrimination is achieved by 
including in the likelihood a variable built, using multivariate analysis tools, from the 
combination (either a linear Fisher discriminant, J-", or a non-linear neural- network, 
NN) of several event-shape quantities. These variables distinguish spherical BB 
events from more jet-like qq events and exploit the different angular correlations in 
the two event categories. B D^*\ decays, which are 12 times more abundant than 
B D^*^K and are expected to show negligible CP- violating effects (r^ ~ 0.01 in 
such decays) , are discriminated by means of the excellent pion and kaon identification 
provided by dE/dx measured in the charged particle tracking devices and by the 
radiation detected in the Cherenkov detector, and are used as control samples. 



that depend on the relative weak phase 7, the magnitude ratio tb = 
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2 Dalitz-plot method: i:)(*)i^(*)±, D Klh+h' 

We reconstruct B^^DK^, D*K^ {D'^D-f and Dn^), and DK*^ (/s:*±->is:°7r±) 
decays, followed by neutral D meson decays to the 3-body self-conjugate final states 
Kgh^h~ {h = 71, K) pj. From an extended maximum likelihood fit to m-^Si and 
J-" (Fig. [T]) we determine the signal and background yields in each channel: we find 
268 B candidates with D K'^K+K^ and 1507 B candidates with D K^tt+tt^. 
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Figure 1: The rngs (a), AE (b), and T (c) distributions for B^-^DK^, D^K'^^n+n- . 
for events in the signal region (mss > 5.272 GeV/c^, \AE\ < 30 MeV, and T > —0.1), 
after all the selection criteria, except the one on the plotted variable, are applied. The 
curves represent the fit projections: signal plus background (solid black lines), qq+BB 
background (dotted red lines), qq + BB + B ^ Dn background (dashed blue lines). 

Following the technique proposed in [2], from a fit to the Dalitz-plot distri- 
bution of the D daughters we determine 2D confidence regions for the variables 
x± = tb cos{6b i 'y) and y± = rBsin(5s±7) (Fig. [2]). In the fit we model the 

and D decay amplitudes to Kgh'^h~ as the coherent sum of a non-resonant part 
and several intermediate two-body decays that proceed through known K^h or h^h^ 
resonances. The model is determined from large (~ 6.2 x 10^) and very pure (~ 99%) 
control samples of D mesons produced in D*^D7t decays [3J. The results for x and 
y are summarized in Table [TJ 
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B± ^ 








x+ 


-0.103 


± 


0.037 


± 0.006 


± 


0.007 


0.147 


± 0.053 


± 


0.017 


± 


0.003 


-0.151 


±0.083 


± 0.029 


± 


0.006 


y+ 


-0.021 


± 


0.048 


± 0.004 


± 


0.009 


-0.032 


± 0.077 


± 


0.008 


± 


0.006 


0.045 


±0.106 


± 0.036 


± 


0.008 


X- 


0.060 


± 


0.039 


± 0.007 


± 


0.006 


-0.104 


± 0.051 


± 


0.019 


± 


0.002 


0.075 


±0.096 


±0.029 


± 


0.007 


y- 


0.062 


± 


0.045 


± 0.004 


± 


0.006 


-0.052 


± 0.063 


± 


0.009 


± 


0.007 


0.127 


±0.095 


±0.027 


± 


0.006 



Table 1: Values of x± and y± measured with the Dalitz-plot analysis of B^^D'^*^ K^*^^ 



From the (a;±, y±) confidence regions we determine, using a frequentist procedure, 
Icr confidence intervals for 7, and 5b (Fig. [3]). We obtain 7 mod 180° = (68 ± 14± 
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Figure 2: Icr and 2a contours in the x±,y± planes for (a) B — )■ DK, (b) B — )• D*K 
and (c) B — )• DK*, for B~ (solid lines) and 5+ (dotted lines) decays. 

4±3)°, where the three uncertainties are respectively the statistical, the experimental 
systematic and the Dalitz- model systematic ones. We find values of around 0.1, 
confirming that interference is low in these channels: r^^ = 0.096±0.029; r^*^ = 
0.133j:[]:J]^^; kr^^'^ = 0.149l[J:[J^^ {k=0.9±0.1 takes into account the K* finite width). 
We also measure the strong phases (modulo 180°): (5g^* = (119t^|^)°; S^'^^ = {-82± 
21)°; 6^^* =(111 ± 32)°. A 3.5cr evidence of direct CP violation is found from the 
distance between (x+,?/+) and (x_,?/_) (0 in absence of CPV) in the three B decay 
channels. 
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Figure 3: 1-confidence level (CL) as a function of 7 (left), tb (center) and 63 (right) 
from the B — ^ Dalitz-plot analysis. 



3 GLW method: DK^*^^, D /(cp) 

We reconstruct B"^ — > DK^ decays, with D mesons decaying to non-CP {D^ — )■ 
K-'K+), CP-even {K^K', Tr+vr") and CP-odd {K^ii^ , K^uj) eigenstates ^. 
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The partial decay rate charge asymmetries Acp± for CP-even and CP-odd D final 
states and the ratios Rcp± of the charged-averaged B meson partial decay rates in CP 
and non-CP decays provide a set of four observables from which the three unknowns 
■y, tb and 6b can be extracted (with an 8-fold discrete ambiguity for the phases) [5]. 

The signal yields, from which the partial decay rates are determined, are obtained 
from maximum likelihood fits to m-Es, Ai? and J-". An example is shown in Fig. |H 
We identify about 500 DK^ decays with CP -even D final states and a similar 

amount of i?^ — t- DK"^ decays with CP-odd D final states. We measure Acp+ = 
0.25 ± 0.06 ± 0.02 and and Acp- = -0.09 ± 0.07 ± 0.02, respectively, where the first 
error is the statistical and the second is the systematic uncertainty. The parameter 
Acp+ is different from zero with a significance of 3.6 standard deviations, constituting 
evidence for direct CP violation. We also measure Rcp+ = 1-18 ± 0.09 ± 0.05 and 
R^p_ = 1.07 ±0.08 ±0.04. 
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Figure 4: AE projections of the fits to the data: (a) B^^Dcp+K~ , (b) 
B^^Dqp+K^ . The curves are the full PDF (solid, blue), and B^Dn (dash-dotted, 
green) stacked on the remaining backgrounds (dotted, purple). We require candidates 
to lie inside a signal-enriched region: 0.2 < F < 1.5, 5.275 < m-Es < 5.285 GeV/c^, 
charged particle from the B passing kaon identification criteria. 

Using a frequentist technique, including statistical and systematic uncertainties, 
we obtain 0.24 < rp < 0.45 (0.06 < tb < 0.51) and, modulo 180°, 11.3° < 7 < 
22.7° or 80.9° < 7 < 99.1° or 157.3° < 7 < 168.7° (7.0° < 7 < 173.0°) at the 
68% (95%) confidence level (Fig. E]). To facilitate the combination of these mea- 
surements with the results of the Dalitz-plot analysis, we exclude the D — )■ K^cj), 
(f) — )■ K^K~ channel from this analysis - thus removing events common to the 
two measurements - and express our results in terms of the variables x± using 
x± = \ [Rcp+i^ T Acp+)-Rcp-il T Acp-)]. We find: x+ = -0.057 ± 0.039 ± 0.015 
and x_ = 0.132 ± 0.042 ± 0.018, in good agreement with the results from the Dalitz- 
plot analysis. 
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Figure 5: 1-CL as a function of 7 mod 180° (left) and rs (right) from the B^DK 
GLW study. 

4 ADS method: D^*^K^, D K^tt^ 

We reconstruct B^^DK^ and D*K^ {D* ^D'-) and /^vr^), followed by D decays to 
both the doubly- Cabibbo-suppressed final state K^n^ and the Cabibbo-allowed 
final state K'tt^, which is used as normalization and control sample [6]. Final states 
with opposite-sign kaons are produced from the interference of the CKM favored 
B decay followed by the doubly Cabibbo-suppressed D decay and the CKM- and 
color- suppressed B decay followed by the Cabibbo-allowed D decay, and the CP 
asymmetries may be potentially very large. On the other hand, their overall branching 
fractions are very small (0(10^^)) and background suppression is crucial. The three 
branching fraction ratios {Rads) between B decays with opposite-sign and same-sign 
kaons and the three charge asymmetries (Aads) in B decays with opposite-sign kaons 
provide six observables that can be used, together with the measurements by c- and 
5-factories of the amplitude ratio r£, and the strong phase difference 6£, between the 
two D decay amplitudes, to determine 7 (with a 4-fold discrete ambiguity) and the 
two sets of tb^Sb [Z|- 

The yields are determined from fits to rngs and NN (Fig. |6]). We see indications 
of signals for the B — > DK and B — )■ D^^oK opposite-sign modes, with significances 
of 2.1(7 and 2.2a, respectively. The measured branching fration ratios are R^ds ~ 
(1.1±0.5±0.2) X 10-2 and i^^^j'Jf = (1.8±0.9±0.4) x lO'^. The CP asymmetries are 
large, A^^g = -0.86±0.47 t^'H and A^j'lf = +0.77±0.35±0.12. We see no evidence 
of opposite-sign B — )■ D*jj^K decays, and measure R^ds = (1-3 ± 1-4 ± 0.8) x 10"^ 
and A^l^s = +0-36 ± 0.94 toU From these results we infer rg^* = O.OQStg;!]!}, 
r^*K± = 0.096toml and 54° < 7 < 83° (Fig. C]). 
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Figure 6: rrj-ES projection of the fit to the data for the — )■ DK^, D — t- K^tt^ 
decays, for samples enriched in signal [NN > 0.94), for (a) and (b) candidates. 
The curves represent the fit projections for signal plus background (solid), the sum 
of all background components (dashed), and the qq background only (dotted). 




Figure 7: 1-CL as a function of 7 (left) and tb (right) from the B—^D^*^K ADS study. 

5 Conclusion 

The full BABAR dataset has been exploited to measure the CKM angle 7 in several 
5± ^ decays using three alternative techniques. A coherent set of results 

on 7 and on the hadronic parameters characterizing the B decay amplitudes has been 
obtained. The central value for 7, around 70°, is consistent with indirect determi- 
nations from the CKM fits. We attained a precision on 7 around 15°, and confirm 
the theoretical expectations of significant suppression (r^ ~ 0.1) of the 6 — )■ m medi- 
ated decay amplitud with respect to the 6 — )■ c one. Finally, two direct CP violation 
evidences at the level of 3.5a have been observed. 
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